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O ■ ABSTRACT 
(N 

' I show that neutral lithium couples strongly to the cosmic microwave background 

(CMB) through its 6708A resonant transition after it recombines at z ~ 500. At 
. observed wavelengths of ^ 335/im, the CMB anisotropics are significantly altered since 

the optical depth for resonant scattering by neutral lithium is substantial, tui ~ 0.5. 
I The scattering would suppress the original anisotropics, but will generate strong new 

\ anisotropics in the CMB temperature and polarization on sub-degree scales {£ ^ 100). 

I Observations at different wavelengths in this spectral regime can probe different thin 

' slices of the early universe. The anisotropy noise contributed by continuum radiation 

f-*) ■ from foreground far-infrared sources could be taken out by subtracting maps at slightly 

different wavelengths. Detection of the above effects can be used to study structure 
at z ^ 500 and to constrain the primordial abundance and recombination history of 
'"q^I lithium. 

o 

^: 

^ ■ Subject headings: Cosmology: theory — cosmic microwave background 



1. Introduction 

The dynamical coupling between the baryons and the Cosmic Microwave Background (CMB) 
during the epoch of cosmological recombination (z ~ 10'^) is commonly assumed to be mediated 
only through Thomson scattering (Peebles & Yu 1970). This coupling becomes ineffective when 
hydrogen recombines, at around the same time when the universe becomes transparent, and hence 
it leaves a clear signature on the temperature anisotropics of the CMB (White, Scott, & Silk 1994, 
and references therein). Forthcoming satellite missions, such as MAPQ in 2001 or Planckg in 2007, 
offer the sensitivity to probe deviations from this simple model at the percent level (Scott 1999; 
Yu, Spergel, & Ostriker 2001). 



^http: / /map. gsfc.nasa.gov 
^http: //astro. estec.esa.nl/Planck 
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In this Letter, I examine the potential significance of the couphng between the cosmic gas and 
the CMB due to resonant Une transitions of neutral atoms. For such transitions, the integral of 
the absorption cross-section over frequency, 

/''('''''■' = ^"(^)> 

is typically many orders of magnitude larger than for Thomson scattering, where /12 is the 
absorption oscillator strength (Rybicki &; Lightman 1979). For example, the Lya transition of 
hydrogen provides an average cross-section which is seven orders of magnitude larger than the 
Thomson value, when averaged over a frequency band as wide as the resonant frequency itself 
[(Av/v) ~ 1]. Since the abundance of atoms in excited states is highly suppressed by large 
Boltzmann factors at the recombination redshift, we consider only transitions from the ground 
state. For helium, these transitions correspond to photon frequencies that are too far on the Wien 
tail of the CMB spectrum and only couple to a negligible fraction of the radiation energy density. 
Hence, we focus our discussion on hydrogen and lithium. 

In §2, we derive the resonant— line opacity and drag force for atoms that move relative to a 
blackbody radiation field. In §3, we consider the cosmological implications of these results for 
hydrogen and lithium. Finally, §4 summarizes the main conclusions of this work. Throughout the 
paper, we adopt the density parameters Qq « 0.3 for matter and i^b = 0.04 for baryons, and a 
Hubble constant Hq 70 km s^^ Mpc^^. The effect of the cosmological constant can be neglected 
at the redshifts of interest here. 



2. Drag Force and Opacity Due to a Resonant Line 

Consider a single atom moving at a non-relativistic peculiar velocity v relative to an isotropic 

blackbody radiation field of temperature T. Because of its motion, the atom sees a dipole 
anisotropy in the radiation temperature, Ta{0) = T[l + {v/c) cos0], where 9 is the angle relative to 
its direction of motion (Peebles 1993). We assume that the atom couples to the radiation through 
a resonant line of frequency uq with (hvo/kT) ^ 1. Hence, the line resonates in the Wien tail of 
the spectrum, for which the radiation energy density per unit frequency (in erg cm~^ Hz~^) is 
given by, 

where we have used the fact that [(hv/kT)(v/c) cos 9] <C 1. Absorption of a resonant photon by 
the atom is followed by isotropic emission in the atom's rest frame, and so the net drag force 
(momentum transfer per unit time) is given by, 

roo rl 

F = - du - dcos9 a{u)u^{9) X cos 9 . (3) 
Jo 2 J-i 
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Since the absorption cross-section is sharply peaked around uq, we may use equation (|^) to get, 

For an ensemble of atoms of density na, which are all moving at the same peculiar velocity v, 
the force per unit volume is given by UgF. If the resonant atoms are embedded within a hydrogen 
plasma of total (free+bound) proton density Hp and are strongly coupled to the plasma through 
binary particle collisions, then the cosmological equation of motion for the gas as a whole reads, 

{rUpUp + mana) + Hv^ = naF, (5) 

or equivalently, 

where H = {a /a) is the Hubble expansion rate, and ta is the characteristic time over which the 
peculiar velocity of the gas is damped due to the drag force on the resonant atoms, 

SnipX^ m.c^ kT (hup] 
- 87r2/i2c/i2 e2 Xa Put/- ^ ' 

Here, Aq = c/uq is the resonant wavelength, rup is the proton mass, and Xa = na/(np + Ag^Ua), 
where Ag, = nig/mp is the atomic weight of the resonant atoms. In the cosmological context, 
T = 2.725 K X (1 + z) is the CMB temperature (Mather et al. 1999). Both T and Xg are functions 
of cosmic time, and so equation @ admits the solution, v = [fo/a(i)] exp {— / dt/ts,{t)}, where 
vq = const and a{t) is the cosmological scale factor. In general, the right-hand-side of equation @ 
should include other terms which source the peculiar velocity field, such as gravitational or 
pressure forces due to density inhomogeneities. The significance of the drag force can be calibrated 
in terms of the Hubble time, H^^{z), through the product i/ta- For the redshifts of interest, 
Hiz) « VnoHo{l + zf/\ 

The assumption of an isotropic radiation background is only valid if the moving fluid element 
is optically-thin to resonant scattering. In the opposite regime of high opacity, the background 
photons are isotropized locally through multiple scatterings and exert a reduced drag force on an 
atom embedded deep inside the moving fluid element. The Sobolev optical depth for resonant 
scatterings in an expanding and nearly uniform universe is given by (Sobolev 1946; Dell'Antonio 
& Rybicki 1993), 

vre^ Aona(z) A2i(52/ffi)Aona(2) , , 

Ta(z = /l2 TTr-r = o rr/ \ ' 

where A21 is the transition probability per unit time for spontaneous emission between the two 
energy levels (in s~^), and {g2/gi) is the ratio between the statistical weights of the excited and 
ground states. The optical depth is dominated by the velocity gradient of the Hubble flow, because 
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gravitationally-induced peculiar velocities perturb this flow only slightly at the early cosmic times 
of interest here. 

For Ta ^ 1, the flux of resonant photons inside the moving fluid element is reduced by a 
factor of ^ Ta (Rybicki & Hummer 1978), and so the drag force in equation is lowered by the 
same factor. This reduction has dramatic consequences for the resonant Lya transition of neutral 
hydrogen, as we show next. 

3. Implications For Cosmological Recombination 
3.1. Hydrogen 

For the Lya transition of hydrogen, Aq = 1216A , A21 = 6.262 x 10® s-^ 

/12 = 0.4162, and (52/51) = 3, yielding an exceedingly large Sobolev optical-depth, 
th = 7.5 X 10®Xh[(1 + 2;)/10^]^/^, where Xh denotes the neutral fraction of hydrogen. The 
thermally-broadened Lya line has an effective cross-section of (a) ~ 3 x 10~^^ cm^ (Peebles 1993) 
and a corresponding mean- free-path of l/nii{a) ~ 10^^ cm, at 2: ~ 10'^. Only on smaller scales 
can Lya viscosity damp velocity gradients effectively. On the much larger scales where the CMB 
anisotropics are measurable, the drag time in equation (^) is lengthened by a factor of ^ th, which 
makes it very much longer than the Hubble time. The situation is similar for deuterium. Balmer 
line transitions of hydrogen have a small optical depth, but the occupation probability of their 
lower level is suppressed by an exceedingly large Boltzmann factor, ~ exp{— 4.35 x 10'^/(1 + z)}. 

The resulting drag force is small, even when including the enhancement in the Lya flux beyond 
the blackbody value, due to the Lya photons released during recombination (~ 1 Lya photon per 
hydrogen atom; see Dell' Antonio & Rybicki 1993 for the distorted spectrum). 

3.2. Lithium 

3.2.1. Drag Force and Opacity 

Lithium has a closed n = 1 shell with two electrons, and one electron outside this shell. The 
transition between the ground state (25) and the first excited state (2P) has Aq = 6708A and 
A21 = 3.69 X 10^ s-\ with /12 = 0.247 and (52/ffi) = 1 for the 2'^S-2'^P^/^ transition, and 
/12 = 0.494, (52/91) = 2 for the 2'^S-2'^Pl^^^ transition (Radzig k Smirnov 1985; Yan & Drake 
1995). We therefore get from equations (0) and (|8|), 

in-4^1 + ^\'^V ^Lil / 7.89 X 10^ 1 

Htui = 4.96 X 10 5- exp < > , (9) 

V 103 y V3.8X 10-10; ^\ l + z j ^' 
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l + z 

103 



) 
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TLiI 



= 2.82 



3.8 X 10-10 



(10) 



Figure 1 depicts the redshift dependence of the drag time in equation (P) and compares it to 
the drag time due to Thomson scattering in the recombining plasma, 




where Xf, = {ue/rip) = ne/(ne + nn) is the ionization fraction of hydrogen. We calculated Xe{z) 
from the updated version of the standard recombination model (Seager, Sasselov, & Scott 2000). 
Despite the low abundance of lithium, its resonant drag force is not much lower than the Thomson 
drag force at z ~ 500, when roughly half of the lithium atoms are expected to have recombined 
(Palla, Galh, & Silk 1995; Stancil, Lepp, & Dalgarno 1996, 1998). 

The net effect of the lithium drag is highly sensitive to the recombination history of lithium 
and hydrogen. Previous discussions (e.g., Stancil et al. 1996, 1998) did not include the delaying 
effect of Lya photons on the recombination history of hydrogen (Peebles 1993), and therefore 
underestimated the free electron fraction and correspondingly the neutral lithium abundance at 
z <^ 10^. Another process with an opposite sign that was omitted, involves the excitation of neutral 
lithium from its ground level by the CMB radiation field. For this effect, it is important to take 
account of the Lyman and Balmer series distortions of the CMB spectrum (Dell'Antonio & Rybicki 
1993), which redshift into resonance at lithium recombination. In fact, the lithium chemistry in 
the early universe might provide a way for inferring the existence of these distortions. Finally, the 
recombination history of '^Be should also be followed since a substantial fraction of the ^Li forms 
through electron-capture in the ^Be produced by Big-Bang nucleosynthesis (T. Walker, private 
communication). Since the half-life of neutral ^Be is only 53 days, this conversion can be regarded 
as instanteneous on the cosmological timescale, as soon as ^Be recombines. Careful calculations of 
the neutral lithium fraction as a function of redshift are required in order to quantify the imprint 
of the lithium drag on the CMB anisotropies (Dalgarno, Loeb, & Stancil 2001). 

The lithium chemistry could also be affected by a large drift velocity between the Li I fluid 
and the hydrogen plasma. To examine this effect, we consider the equation of motion for the Li I 
fluid (Burgers 1969), 



Here Wui = S^^MkiLiW^kiLi is the collision rate of a Li atom with the gas, where Wk:Lii = n]^{o'k:Ui'v) , 
and ^k;Li = ['^T'ki^Li/i^pinT'k + "iLi)]- The index k runs over all other particle species in the gas 
which are assumed to move at a common bulk velocity v because they are not subject to the 
resonant drag force. The geometric cross-section of Li I atoms, ciLii ~ 10"^^ cni'^, couples them to 
the rest of the gas on a time scale ~ 10'' s, much shorter than the Hubble time, ~ 3 x 10^^ g 




(12) 
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at z ~ 10^. Hence, we may ignore the terms on the left-hand-side in equation (12), as well as the 
perturbed gravitational force which is of the same order. The remaining terms give the solution, 

V ~1 + Y ^^"^^ 

where Y = WuituiXui/Au. Substituting Au = 7, we find that Wui ~ 10~'^[(1 + z)/10^]^-^ s'^, 
and that Y ^ 1 at z ^ 400, after lithium recombination. At higher redshifts, the assumption that 
vui = V seems not to apply and one may wonder whether the drag time in equation (^) should 
be multiplied by a factor of (v/vui) = (1 + ^)/^ > 1) since the drag force acts on VLii and not 
on V. In this context it is crucial to realize that the Li I nuclei do not maintain their identity 
as neutral atoms; in fact, they continuously form and get ionized on a timescale much shorter 
than a Hubble time. The cross-section for Coulomb collisions of Li+ is far greater than that of Li 
I, au+ = 2y/l^e^\nK/{kTf = 2.7 x 10"^^ [(1 + zj/lO^]-^ cm^, where In A « 20 is the Coulomb 
logarithm. The fraction of time spent by a Li nucleus in an ionized state is the ionization fraction, 
/+ = A'Li+/(XLii + ^Li+)i ^iid so (TLii should be replaced by [fj^ai^i+ + (1 — /+)(TLii] and the drag 
force term should be multiplied by (1 — /+) in equation (^). We consequently find that the 
substantial ionization fraction of the lithium fluid (Standi et al. 1996, 1998) diminishes its drift 
velocity relative to the gas at all relevant redshifts, z <^ 10^. Hence, non-thermal relative velocities 
may be ignored in the Li I chemistry. 

At z ^ 10^, the mean-free-path for resonant scattering by lithium exceeds the relevant 
spatial scales for the CMB anisotropies (10^^-10^^ cm), even if the lithium abundance is taken 
to be somewhat higher than the value, Xsbbn = 3.8 x 10~^°, predicted by the latest deuterium 
measurements and Standard Big-Bang Nucleosynthesis (Buries, Nollett, & Turner 2001; see also 
Walker et al. 1991; Smith, Kawano, & Malaney 1993). Note that values as high as ^Lii ~ 10~^ 
were suggested by models of inhomogeneous Big-Bang nucleosynthesis (Applegate & Hogan 1985; 
Sale & Mathews 1986; Mathews et al. 1990), and would result in Hti^n <^ 10 for velocity gradients 
on sufliciently small scales with rm <^ 1. 



3.2.2. Effects on CMB Anisotropies 

If a substantial fraction of lithium recombines by z ~ 500 (as predicted by Palla et al. 1995 
or Stancil et al. 1996, 1998), then the flux of the original anisotropies will be suppressed by the 
absorption factor, exp(— TLii), at observed wavelengths below (6708A x500) = 335/um, for which 
''"Lii ~ 0.6(XLii/A'sBBN)- The resonant scattering effect is different from that caused by Thomson 
scattering at the epoch of reionization (Hu & White 1997), in that it has a larger optical depth and 
it occurs within a much thinner shell of gas (for a sufficiently narrow band of observed wavelengths) . 
Consequently, it should induce new first-order anisotropies [multiplied by (1 — exp(— tlu)] due 
to the coherence of the velocity field in the thin scattering shell. The dominant contribution 
to the new anisotropies would come from the Doppler effect. In contrast to reionization, there 
should be no severe cancellations of the Doppler effects from line-of-sight velocity fluctuations at 
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large wavenumbers, because the last scattering surface (or "visibility function") for the resonance 
is extremely thin. The increase in amplitude of sub-horizon fluctuations between z = 10^ and 
z = 500 would enhance the new anisotropics relative to the original ones and introduce a distortion 
of order unity in the original anisotropy spectrum on these scales. The anisotropy spectrum 
would be modified at multipole moments i ^ 100, around and below the angular scale of the first 
acoustic peak (which is a factor of ~ -v/6 = 2.4 smaller than the horizon scale at z ~ 500), where 
the original anisotropy amplitude reaches its maximum (i ~ 220). Resonant scattering would also 
result in enhanced polarization anisotropics on the same angular scales (Chandrasekhar 1960). 
Since the scattering is done by lithium, the anisotropics will reflect the peculiar velocity of the 
lithium fluid which, as shown above, is not expected to deviate significantly from the velocity of 
the whole. 

Since the temperature fluctuations are in the Wien tail of the CMB spectrum, they translate 
to brightness fluctuations (in erg s~^cm~^sr~^Hz~^) of a much larger contrast, 



where z is the redshift being probed and we have substituted By{T) oc exp^—hv/kT). Hence the 
first acoustic peak would correspond to brightness fluctuations of up to ~ 5 x 10~^. 

At a given observed wavelength, the redshift thickness of the scattering shell is only determined 
by two contributions: (i) the possibility that a photon at this wavelength was either absorbed by 
the 2'^S-2'^P^^2 transition with Ao = 6707.76A or at a slightly lower redshift by the 2'^S-2'^P^/^ 
transition with Aq = 6707.9lA (Radzig & Smirnov 1985), yielding [Az/{1 + z)] = 2.2 x 10~^; and 
(ii) the thermal width of the line^, [A2:/(l + z)] = (AA/A)th ~ 10^^. In practice, the minimum 
thickness would be limited by the band-width of the detector. The maximum band-width and 
minimum angular resolution should be tuned to match the expected coherence length of the 
peculiar velocity field at z ~ 500, which corresponds to (AA/A) ~ 0.2 and A9 ~ 6' or equivalently 



Observations at different far-infrared (FIR) wavelengths would probe different thin slices of 
structure in the early universe. The anisotropy pattern would vary gradually as a function of 
wavelength and "3D tomography" is in principle possible, by which one may probe correlations 
in the structure along the line-of-sight. The feasibility of this measurement might be severely 
compromised by the anisotropics of the FIR foreground. Unfortunately, 350/xm is just the 
wavelength where the CMB intensity blends into the extragalactic FIR foreground (see Fig. 1 
in Scott et al. 2001). However, the contribution of discrete sources (see Fig. 6 in Knox et al. 
2001) could be separated out through observations at shorter wavelengths with a higher angular 
resolution. Furthermore, continuum foreground emission would result in similar anisotropics for all 



^Since tlii < 1, the photon does not need to redshift away from the hne center by more than the hne width in 
order to escape to infinity. 




(14) 



i ~ 10^. 
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wavelengths within a band width (AA/A) 1, while the CMB anistropies would vary considerably 
across such a band since its width corresponds to Az/(1 + z) 1. Subtraction of anisotropy 
maps at slightly different wavelengths within such a band can therefore be used to isolate out the 
CMB component (as long as these wavelengths do not overlap with strong emission lines from 
the foreground sources). For i ^ 100, the anisotropics should have the same power-spectrum as 
measured by MAP or Planck at much longer wavelengths. The predictable correlations between 
the temperature and polarization patterns of the CMB and their general statistical properties 
could be used to further enhance the signal-to-noise ratio (Zaldarriaga & Loeb 2001). Obviously, 
these difficult observations need to be made through holes in the Galactic dust emission at high 
Galactic latitudes. 

4. Conclusions 

The 6708A resonant transition of neutral lithium provides a previously unexplored coupling 
between the baryons and the CMB after cosmological recombination. The lithium opacity is 
substantial at observed wavelengths ^ 335[(1 + zui)/500]fiin, where zui is the redshift at which 
^ 50% of the lithium recombines. At these wavelengths, the CMB anisotropics on sub-degree 
angular scales would be significantly altered, due to the finite optical depth for resonant scattering, 
TLii « 0.5[(1 + 2:Lii)/500]'^/^. This scattering would generate new temperature and polarization 
anisotropies at mutipole moments i ^ 100. The detection of these anisotropies depends on the 
prospects for cleaning the expected anisotropies of the FIR foreground (Knox et al. 2001), and may 
be feasible through the subtraction of anisotropy maps for pairs of slightly different wavelengths. 

The relevant wavelength range overlaps with the highest frequency channel of the Planck 
mission (352^m) and with the proposed baloon-borne Explorer of Diffuse Galactic Emissions^ 
(EDGE) which will survey 1% of the sky in 10 wavelength bands between 230-2000;um with 
a resolution ranging from 6' to 14' (see Table 1 in Knox et al. 2001). However, a new 
instrument design with multiple narrow bands (AA/A ^ 0.1) at various wavelengths in the range 
A = 250-350/xm and with angular resolution ~ 6', is necessary in order to optimize the detection 
of the lithium signature on the CMB anisotropies. 

The effects discussed in this Letter depend critically on the recombination history of lithium, 
and emphasize the need to explore the related lithium chemistry in more detail in the future. 

The author thanks Alex Dalgarno, Daniel Eisenstein, George Rybicki, and Matias Zaldarriaga 
for useful discussions. This work was supported in part by NASA grants NAG 5-7039, 5-7768, and 
by NSF grants AST-9900877, AST-0071019. 
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Fig. 1. — Redshift evolution of the dynamical drag time due to the Li 6708 A resonance with the 
microwave background. The drag time shown by the solid line should be divided by the neutral 
fraction of the primordial lithium abundance (^sbbn = 3.8 x 10~^°, see Buries et al. 2001), as 
a function of redshift. The dotted line shows the result for one of the recombination histories of 
lithium in the literature (Palla et al. 1995). For comparison, the Thomson drag time is plotted as 
the dashed line, taking account of the full recombination history of hydrogen. In all cases, the drag 
times are normalized by the Hubble time, H~^{z). 



